Naval  Research  Laboratory 

Washington,  DC  20375-5320 


NRL/MR/6750— 08-9122 


Preliminary  Global  Radiation  Belt 
Formation  and  Prediction  Model 


Michael  J.  Keskinen 

Charged  Particle  Physics  Branch 
Plasma  Physics  Division 


June  25,  2008 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway, 
Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 
25-06-2008 


2.  REPORT  TYPE 

Memorandum  Report 


3.  DATES  COVERED  (From  -  To) 


4.  TITLE  AND  SUBTITLE 


Preliminary  Global  Radiation  Belt  Formation  and  Prediction  Model 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 


Michael  J.  Keskinen 


5d.  PROJECT  NUMBER 

67-9387-F7 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 
4555  Overlook  Avenue,  SW 
Washington,  DC  20375-5320 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


NRL/MR/6750— 08-9 1 22 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 
One  Liberty  Center 
875  North  Randolph  St. 

Arlington,  VA  22203-1995 


10.  SPONSOR  /  MONITOR’S  ACRONYM(S) 

ONR 


11.  SPONSOR  /  MONITOR’S  REPORT 
NUMBER(S) 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

A  computational  model  for  the  time-dependent  dynamics  of  the  electron  and  ion  terrestrial  radiation  belts  has  been  developed.  The  model 
solves  the  fundamental  bounce-averaged  electron  and  ion  modified  Boltzmann  Fokker-Planck  equations  for  energy,  pitch  angle,  and  L-shell.  The 
code  has  been  applied  to  broadband  whistler  turbulence  and  found  to  be  consistent  with  empirical  models  of  radiation  belt  dynamics. 


15.  SUBJECT  TERMS 

Radiation  belt 
Space  weather 


Satellite  damage 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

18.  NUMBER 

19a.  NAME  OF  RESPONSIBLE  PERSON 

OF  ABSTRACT 

OF  PAGES 

Michael  J.  Keskinen 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

UL 

26 

19b.  TELEPHONE  NUMBER  (include  area 

Unclassified 

Unclassified 

Unclassified 

code) 

(202)  767-3215 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


1.  Introduction 


The  earth’s  radiation  belts  are  composed  of  approximately  tens  of  keV  to  several 
MeV  electrons  and  hundreds  to  several  hundred  MeV  ions  which  are  trapped  in  the 
magnetosphere  from  approximately  1.2  <  L  <  8.  Depending  on  energy,  the  electron 
population  typically  reside  in  two  distinct  inner  and  outer  belt  regions  separated 
by  the  slot  region  (1.8  <  L  <  3.2). 

The  terrestrial  radiation  belts  can  have  a  significant  effects  on  space  systems 
and  space  technology.  Moderate  energy  (10-100  keV)  electrons  can  lead  to  sur¬ 
face  charging  in  space  systems  while  relativistic  MeV  electrons  can  cause  deep- 
dielectric  charging  in  spececraft  materials.  During  geomagnetic  substorms,  large 
fluctuations  in  electron  differential  flux  can  occur  on  relatively  short  time  scales. 
A  time-dependent  predictive  model  for  the  earth’s  radiation  belts  will  have  both 
scientific  and  practical  importance. 

Trapped  charged  particles  in  the  earth’s  radiation  belts  have  multiple-periodic 
motions  on  three  distinct  timescales.  To  describe  this  motion  [ Northrop ,  1963;  Roed- 
erer,  1970]  three  adiabatic  invariants,  i.e,  J i ,  J2,  and  J3,  have  been  developed  with 
Ji  =  7rcp2sin2a/2qB,  J2  =  / Pyds,  J3  =  (q/c)<l>.  Here,  a  is  the  pitch  angle,  p  is  the 
total  momentum,  py  is  the  component  of  momentum  parallel  to  the  geomagnetic 
field,  and  $  is  the  magnetic  flux  enclosed  through  the  drift  shell.  An  alternative  set 
of  invariants  are  M,  J  =  J2,  and  $  with  M  =  p2sin2o;/2moB  with  m0  the  rest  mass. 
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When  turbulence,  e.g.,  from  wave  particle  interactions,  produce  fluctuations  in  the 
J’s  that  are  small,  random,  and  high  frequency,  the  changes  in  f  can  be  modeled  as 
a  stochastic  diffusive  process  and  can  be  described  by  a  Fokker-Planck  formalism. 
In  this  formalism,  the  description  of  radiation  belt  dynamics  involves  a  balance  be¬ 
tween  diffusive  transport,  convection,  sources,  and  losses.  During  magnetic  storms, 
external  sources  of  both  electrons  and  ions  are  important  and  need  to  be  incor¬ 
porated.  Key  sources  of  phase  space  transport  include  radial  (L)  diffusion  due  to 
fluctuations  in  the  convection  electric  field  and  magnetospheric  magnetic  field,  pitch 
angle  and  energy  diffusion,  e.g.,  due  to  whistler  waves,  pitch  angle  diffusion  and 
energy  losses  from  collisional  processes  involving  cold  plasmasphcric  plasma  and 
neutral  atmospheric  species,  and  charge  exchange  losses.  Precipitation  losses  can 
occur  when  the  pitch  angle  is  inside  the  loss  cone. 

Several  models  have  been  developed  to  simulate  global  radiation  belt  dynamics. 
Two  different  approaches  have  been  taken,  i.e.,  the  test-particle  and  kinetic  formu¬ 
lations.  The  test  particle  approach  involves  the  tracing  of  individual  particle  orbits 
in  the  combined  electric  and  magnetic  fields  derived  from  an  MHD  model  of  the 
magnetosphere.  Using  the  test  particle  method,  guiding  center  simulations  [Li  et 
al.„  1993;  Hudson  et  al,  1996]  of  the  sudden  storm  commencement  induced  by  the 
interplanetary  shock  of  March  24,  1991  were  performed  and  analyzed.  These  studies 
suggested  the  formation  of  a  new  electron  belt  at  L=2.5.  In  addition,  test  particle 
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simulations  of  the  outer  belt  electrons  [ Ukhorskiy  et  al,  2006]  were  performed  to 
model  the  magnetic  storm  of  September  7,  2002. 

Kinetic  models  have  also  been  used  to  model  electron  radiation  belt  dynam¬ 
ics.  This  approach  solves  the  modified  Boltzmann  Fokker-Planck  equation  with 
wave-induced  and  collisional  pitch  angle  and  energy  diffusion  coupled  with  parti¬ 
cle  sources  and  losses.  Radial  diffusion,  combined  with  acceleration  and  loss  due 
to  whistler  mode  VLF  chorus,  indicated  that  [ Varotsou  et  al.,  2005;  Horne  et  al, 
2006]  wave  acceleration  by  whistler  waves  is  an  important  mechanism  for  the  outer 
radiation  belt  electrons.  Beutier  and  Boscher  [1995]  and  Bourdarie  et  al.  [1997] 
solved  for  a  phase  space  three-dimensional  model  of  the  electron  radiation  belt. 
They  found  that  cosmic  ray  albedo  neutron  decay  internal  sources  are  not  sufficient 
to  account  for  the  typical  quiet  time  electron  fluxes.  A  convection- diffusion  model 
[Fok  et  al,  2001;  Zheng  et  al,  2003;  Fok  et  al,  2008]  has  been  developed  to  simu¬ 
late  electron  radiation  belt  dynamics.  The  model  was  used  to  simulate  a  substorm 
injection  during  a  dipolarization  of  the  magnetic  held  during  a  substorm.  Experi¬ 
mentally  observable  features  such  as  drift  echoes  were  successfully  reproduced. 

In  this  paper  we  develop  a  first  principles  time-dependent  simulation  model  for 
the  electron  and  ion  radiation  belt  dynamics.  This  is  accomplished  by  solving  the 
convection- diffusion  model  of  radiation  belt  dynamics  as  derived  from  a  modified 
Boltzmann  Fokker-Planck  model.  The  model  includes  both  pitch  angle  and  energy 
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diffusion  from  both  wave-particle  and  collisional  effects,  energy  loss  from  plasma- 
spheric  plasma  and  neutral  atmospheric  constituents,  radial  diffusion,  and  charge 
exchange  losses.  We  find  that  the  general  characteristics  of  the  electron  fluxes  pre¬ 
dicted  by  the  model  are  consistent  with  empirical  models  derived  from  satellite 
observations. 

The  outline  of  the  paper  is  as  follows.  In  Section  2  we  present  the  time- 
dependent  radiation  belt  simulation  model.  In  section  3  we  apply  the  model  to 
whistler  wave  turbulence  and  compare  the  model  output  with  empirical  satellite 
models.  Finally  in  section  4  we  summarize  our  results. 


2.  Model 

Using  kinetic  theory,  for  small,  randomly  phased  fluctuations  in  the  three  adia¬ 
batic  invariants,  J i ,  J2,  and  J3,  the  time  dependence  of  the  electron  and  ion  phase 
space  density  f*.,  k=e,i,  can  be  written  in  the  form  of  a  Fokker-Planck  equation: 


dfk 

dtt 


1 


d_ 

dJi 


3  3  d 

g.a  yy  jprD.i 

i=  1  j= 1  UJi 


^  +  sk-  Lk 


(1) 


where  G  models  energy  loss,  D  describes  diffusive  processes,  and  S*,  and  Lk  represent 
explicit  electron  and  ion  sources  and  losses  within  the  phase  space  volume.  In 
analogy  with  a  Fokker-Planck  model,  G  represents  a  nonstochastic  friction  and  D 
a  stochastic  diffusion.  The  evaluation  of  the  diffusion  tensor  elements  is  of  crucial 


importance  in  the  use  of  Eq.  (1)  for  the  modeling  of  the  radiation  belts.  For  pitch 
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angle  and  energy  diffusion,  both  wave-particle  and  low  altitude  collisional  effects 
will  be  important  while  for  spatial  diffusion,  e.g.,  radial  diffusion,  slow  time  scale 
electric  and  magnetic  fluctuations  will  be  needed. 

The  standard  approach  for  studying  the  effects  of  experimentally  observed  broad¬ 
band  wave  turbulence  on  the  adiabatic  invariants  M,J,<f>  is  to  model  the  effects  of 
wave-particle  interactions  on  pitch  angle  and  energy  diffusion  using  quasilinear  the¬ 
ory  [Kennel  and  Petschek,  1966].  This  theory  was  applied  to  electron  pitch  angle 
scattering  by  whistler  turbulence  by  Lyons  et  al.  [1972]  and  Lyons  and  Thorne 
[1973].  These  studies  established  cyclotron- resonant  whistlers  as  a  key  mechanism 
for  electron  loss  in  the  slot  region  and  for  maintaining  the  steady  state  electron 
distribution  in  combination  with  radial  diffusion  and  Coulomb  collisional  effects  at 
low  altitudes. 


In  this  study  we  include  only  radiation  belt  electron  dynamics.  The  ion  dynam¬ 
ics  is  also  governed  by  a  similar  Fokker-Planck  equation  with  different  pitch  angle 
diffusion  coefficients  along  with  different  source  and  loss  terms. 

It  is  convenient  to  transform  Eq.  (1)  from  (M,  J,  <f>)  space  to  (E,  x,  L)  space 
where  E  is  the  energy,  x=cos  «o  with  a0  the  equatorial  pitch  angle,  and  L  denotes 
L  shell.  Since  we  will  be  interested  in  time  scales  much  longer  than  the  electron 


bounce  time,  we  bounce-average  Eq.(l)  for  the  electrons: 


d(f)  dLd(f)  dcj>d(f)  1  a  d_E 
dt  dtdL  dt  <90  ^PdE  dt)[J) 


1  d 
xT  ( y )  dx 


xT  (: y )  (Dxx) 


m 

dx 
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1  d  T(  wn  tn  \^J1 

~xT(y)dxX  ^  xE)  dE  +7 pdE1P^  Ex)  dx 

1  a_wn,  M  ■  r2*  1  ,  »</)  .  o  , 


In  Eq.  (2),  the  angle  brackets  denote  bounce  averaging,  i.e.,  applying  the  op¬ 
erator  tb  /  ds/vy  with  tb  the  bounce  period  and  the  coordinate  s  denotes  distance 
along  the  geomagnetic  held.  In  Eq.  (2),  the  second  and  third  terms  on  the  left  hand 
side  represent  radial  and  azimuthal  convection,  the  fourth  term  represents  electron 
energy  loss  due  to  collisional  effects  in  the  plasmasphere  and  atmosphere,  the  first 
term  on  the  right  hand  side  represents  pitch  angle  diffusion  from  both  collisional 
effects  and  wave-particle  interactions,  the  second  and  third  terms  represent  cross 
energy  and  pitch  angle  diffusion,  the  fourth  term  on  the  right  hand  side  represents 
energy  diffusion  from  both  collisional  and  wave-particle  effects,  and  the  fifth  term 
gives  radial  L-diffusion.  In  addition,  we  have  defined  y  =  (1  —  x2)1/2  =  sin  «0 ,  and 
T(y)  =  1.38-0.31(y  +  y1/2). 

The  energy  loss  rate  for  the  radiation  belt  electrons  comes  primarily  from  in¬ 
elastic  collisions  with  neutral  atmospheric  atoms  and  elastic  collisions  with  plas- 
maspheric  and  ionospheric  electrons.  The  bounce- averaged  energy  loss  rate  for 
electrons  can  be  written  [Walt  and  MacDonald,  1964]: 


,dE  47 re4  E  +  E0  rx™  sin2 a  cos7 A 

(-7T  = - - r  /  d A - /  (ne,  n,  E) 

dt  mcl  (y)  (E  +  2E0)^  cosasmza 0 
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f  E(-g-  +  2V  1  e 

f  (ne,  n,  E)  =  nJnrjj-,1  +  ^  In  V  0 - }- - -  (4) 

;  J-i  Z  ±j  +  ±j0 

V  / 

In  Eq.  (3),  a  denotes  the  local  pitch  angle,  A  is  the  magnetic  latitude,  Am  is  the 
magnetic  mirror  latitude,  E0  =  me2,  I;  is  the  ionization  energy,  r/D  =  h/mrvAD, 
mr  is  the  reduced  mass,  ne  is  the  electron  density,  n,  is  the  neutral  density  of 
species  i,  and  Z  is  the  atomic  number.  The  electron  and  ion  densities  H+,  He+,  0+, 
O2",  NO+  are  found  using  the  International  Reference  Ionosphere  [ Bilitza ,  2000] 
combined  with  a  plasmaspheric  model  [ Chiu  et  al,  1979].  The  neutral  species  are 
found  using  the  MSIS90  model [Hedin,  1991]. 

The  bounce-averaged  pitch  angle  diffusion  coefficient  (D**)  in  Eq.  (2)  is  com¬ 
posed  of  a  collisional  part  and  wave-particle  interaction  part  and  can  be  written 
Dxx  =  D]]x  +  Dxx.  Similarly,  the  bounce-averaged  energy  diffusion  coefficient  (Dee) 
is  composed  of  a  collisional  part  and  wave-particle  interaction  part  and  can  be 
written  DEE  =  D|E  +  DgE.  The  cross  tensor  DxE  and  DEx  elements  are  produced 
primarily  by  wave-particle  interactions.  The  bounce-averaged  collisional  pitch  angle 
scattering  diffusion  coefficient  can  be  written  [  Walt  and  MacDonald,  1964]: 


(Dxx)  = 


27re4c  (E  +  E0)  rxm  cosa  cos7 X 
El  (E  +  2E0) 3  T  (y)  Jo  cos2a0 


g(ne,n,E) 


g  (ne,  n ,  E)  =  nplngDl  +  nelngDl  +  ^  Z'-/rillnril  1 


(6) 


Tji  =  Z1//3(l  — /32)1/2/137/9,  and  f3  =  v/c.  The  collisional  energy  diffusion  coefficient 
is  typically  much  smaller  than  the  collisional  pitch  angle  diffusion  coefficient  and 
will  not  be  considered  further. 

In  order  to  estimate  the  wave-particle  contribution  to  pitch  angle  and  energy 
diffusion,  a  specific  wave-particle  interaction  must  be  specified.  Wave-particle  inter¬ 
actions  play  a  fundamental  role  in  radiation  belt  electron  dynamics.  Gyroresonant 
interactions  are  produced  by  the  local  Doppler-shifted  cyclotron  resonance  condi¬ 
tion,  u  —  sDe/7— kvcoso;  =  0.  Here  s=l  denotes  right  hand  R-mode  waves  and  s=-l 
the  left  hand  L-mode  waves  and  7  is  the  Lorentz  factor.  Gyroresonant  interactions 
can  occur  with  both  high  and  low  frequency  waves  in  the  range  of  approximately 
0.1f2o+  <  uj  <  0.8De  where  Do+  is  the  oxygen  ion  gyrofrequency  and  De  is  the  elec¬ 
tron  gyrofrequency.  Waves  in  this  frequency  range  include  whistler  mode  chorus, 
plasmaspheric  hiss,  and  electromagnetic  ion  cyclotron  (EMIC)  waves. 

In  this  study  we  consider  only  R-mode  whistler  turbulence.  The  evolution  of  f  in 
the  presence  of  whistler  wave  turbulence  can  be  found  from  the  standard  quasilinear 
theory  [ Lyons  et  al,  1974]: 


df  1  d  .  df  : 

■777  :  "o  DaQsifiQL—  b  :  x 

at  srna  aa  oa  svna  001 


d  n  ■  df 
— uavsvna— 

op 

d 


1  d  2  df  1  d  2  df 
p2  dpP  pa  da  p2  dpP  pp  dp 


(7) 


For  parallel  propagating  R-mode  whistler  waves  the  pitch  angle  diffusion  coef- 
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ficient  can  be  written  [ Summers ,  2005]: 


D  =  ^ 


El 


E 


( 1  -  K~vCOSa )' 


2  Bl  (E  +  E0f  i  iy  cosa  -  vg) 


■I  {hi 


(8) 


The  cross-diffusion  coefficient  can  be  written: 


n  _  71  Q?P  E0sina  v-  fe)  (X  ~~  MC0S“)  r/,„ 

-^“P  O  D2  /  771  ,  77  \2  2—d  / _  ..  \  J  Vfi'* 


2  (£  +  E{ 


0) 


(■ vcosa  —  vg) 


(9) 


The  energy  diffusion  coefficient  can  be  written: 


p  ^eP2  El  sin2  a  v  (m) 

^PP  ~  r,  D2  /n  n  L 


2  (E  +  £0)2  i  (^cosa  -  Vg) ' 


(10) 


In  Eq.  (8)-(10)  the  a 4  and  k,:  satisfy  both  the  cyclotron  resonance  condition  with 
k  =  kl! 

fle  1  i 

uji  —  hi  v  cosa - =  0  (11) 

7 

and  the  cold  plasma  linear  dispersion  relation: 


c2k2 


cot 


h  =  1 


u, 


pe 


(idi  (j^i  “I-  ^i) 


(12) 


The  whistler  wave  spectrum  I(k)  in  Eq.  (8)-(10)  is  defined  by: 


B?  rK  2 

s i=L  I(k)dk 


rk  2 

I  k\ 


(13) 


where  the  total  magnetic  field  B  =  B0  +  Bi  is  composed  of  a  background  part  B0 


and  a  small  fluctuating  part  B| . 
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Including  whistler  wave  convection,  the  quasilinear  equation  for  the  whistler 
wave  spectrum  I  can  be  written: 


9  I  Ouj  d  I  dijj  d  I 
dt  ^  9k  dr  dr  9k 


2y  (k)  I  (k) 


(14) 


where  vg  =  dco/dk  is  the  whistler  group  velocity  and  the  the  growth  rate  7  given 
by  [Kennel  and  Petschek,  1966]: 
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2vr2fL 


vrV  (■ vR ) 


A  (vR) 


(15) 


with 


77  (vR)  =  27 tvr  dv±v±f  (y±,  uy  = 


(16) 


and 


\  1  \  1  f°°  J  V1  (  9/  9/ 


V\\=VR 


(17) 


Here  vr  is  the  resonant  velocity  determined  from  Eq.  (11). 

Transforming  from  (v_l,V||)  coordinates  to  energy  and  pitch  angle  coordinates 
(E,  a)  the  growth  rate  can  be  written: 

1 


7  =  2^(l-^7^V^C(S«) 


a 


rnP 


Eft 


Ae  (Er) 


(18) 


with 


C  (Ei;!  =  j  da  tana  f  (a,  ER) 


(19) 


and 


Ae  (Er) 


1 _ 2ErE0  r^tan2adf  ( a,ER ) 

C  ( Er )  Er  (Er  +  2 Eq)  J  a  cos2a  da 


(20) 
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where  Er  =  (l/2)mev^. 

Consistent  with  experimental  observations,  we  consider  a  whistler  wave  spec¬ 


trum  of  the  form: 


B 2  1  _t  u-um. ' 


7T 2  p  ( LO-m  ^2  \  p  ( ^  2  ^  m 

P  =  TT  erf  - 7 -  +  er'f  - 7 - 


Here  Bx  is  the  rms  wave  magnetic  held,  loi  is  the  lower  frequency  limit  to  the 
wave  spectrum,  U2  is  the  upper  frequency  limit  to  the  wave  spectrum,  ujm  is  the 
center  frequency,  5u  is  the  frequency  bandwidth  and  erf  is  the  error  function.  To 
match  typical  experimental  observations,  uj\  =  0.05h2e,  ca2  =  0.65f2e,  u>m  =  0.35h2e, 
5uj  =  0.15h2e,  Bi  =  lOOpT,  and  c oe/Vte  =  2.5.  We  have  I(cn)  =  (27t/vs)I (k).  Using 


the  form  for  the  wave  spectrum  in  Eq.  (21)-(22)  we  can  write  Eq.  (8)-(10): 


p  El 

J-S  rvrv  ~ 


BA2!1 


2  p  (E  +  E0)  t  \B0J  5uj  (v  cosa  -  vg) 


k  COSaj  Vg  /Uj-UJm  \2 


_  g  \ 


The  cross-diffusion  coefficient  can  be  written: 


D 


n<Pl2ep  El sina  ^  /-BA2  (lEcosa)  vg 


2  (3p  (E  +  Eq)  ;  v60/  Su  ( vcosa  -  vg) 


The  energy  diffusion  coefficient  can  be  written: 


Dpp 


ntijp2  El  sin2  a  y,  /HA2  {jEcosa)  v9  „-A’AmV 
2  f32p  (E  +  E0)2  i  \B0)  8u  ( vcosa  —  v0)6 
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For  inclusion  into  Eq.  (2),  the  bounce- averaged  pitch  angle  diffusion  coefficient 
can  then  be  written: 

=  A  (26) 

1  [y )  Jo  cosza0 

In  addition,  the  bounce-averaged  energy  diffusion  coefficient  can  be  written: 

=  1  (27) 
w  1  ( y )  Jo  svnz cto  cosck 

The  bounce-averaged  cross  diffusion  coefficient  can  be  written: 

1  C»ii!WA 

F  i  (y)  Jo  svnao  cosao 

Including  whistler  wave  convection  but  ignoring  inhomogeneous  plasma  effects, 
the  quasilincar  equation  for  the  bounce- averaged  whistler  spectrum  I  in  Eq.  (14) 


can  be  written: 


1,(1)  =  27(*0  (I) 


where  =  r”1ln(l/R),  rg  =  LRe/v9,  and  R  is  the  ionospheric  reflection  coefficient. 
For  perfect  reflection,  R=1  while  for  strong  ionospheric  absorption,  R  <<  1.  Here 
7 g  models  wave  convection  effects  out  of  the  magnetic  flux  tube. 


In  Eq.  (2)  we  have: 


c*(E  +  2E„) 

-  (E  +  E0f  {°rr) 


c2E 

(DaE)  =  (Dap) 

+  ±JQ 


(31) 
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For  the  convection  terms  in  Eq.  (2)  we  assume  that  we  can  write  the  total  drift 
velocity  as 


V,4(.tMA)BxVS  +  ^  (32) 

with  E  =  Ecor  +  ~EiConv  with  Ecor  and  E conv  the  corotation  and  convection  electric 
fields,  respectively.  We  can  then  write: 

d<t>  —Qnp2c  D{y )  ,(df\  /oox 

dt  mL2R2eeBy2T(y )  cor  \dt )  conv  [  } 


and 


dL 

dt 


(34) 


with  D(y)  =  0.45  -  0.19y  -  0.106ys/4. 

3.  Results 

For  this  preliminary  study,  we  solve  Eq.  (2)  numerically  including  only  en¬ 
ergy  loss,  collisional  pitch  angle  diffusion,  and  whistler  wave  pitch  angle  diffusion 
together  with  a  source  and  loss  term.  Radial  and  azimuthal  convection  as  well 
as  radial  diffusion  effects  are  not  included.  The  total  pitch  angle  diffusion  co¬ 
efficient  is  the  sum  of  Eq.(5)  and  Eq.(26).  For  the  source  term  we  use  Se  = 
(me/2E0)3/2(l/n0)(dne/dt)sin2a:0exp(— E/E0).  For  the  loss  term  we  take  Le  = 
4g(a;o)(f )/tb  with  g(cKo)  =  0  for  pitch  angles  outside  the  loss  cone,  i.e.,  ao  >  a l 
and  g(ao)  =  1  f°r  pitch  angles  inside  the  loss  cone,  i.e.,  ao  <  ckl-  We  treat 
the  range  of  L-shcll,  pitch  angle,  and  energy  of  1  <  L  <  7,  3°  <  a  <  90°,  and 
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0.1  MeV  <  E  <  5MeV. 

Fig.  1  shows  the  bounce-averaged  pitch  angle  diffusion  coefficient  Eq.  (26)  as 
a  function  of  equatorial  pitch  angle  for  several  energies.  Here  we  have  used  the 
whistler  wave  spectrum  Eq.  (21)  for  parallel  propagating  whistler  turbulence  with 
mean  wave  amplitude  of  B=100  pT  situated  at  L=4.5.  The  pitch-angle  diffusion 
coefficient  for  the  MeV  electrons  is  smaller  compared  to  the  smaller  energies. 

Fig.  2  gives  an  example  of  the  omnidirectional  differential  electron  flux  at  a 
fixed  energy  of  0.5  MeV  in  the  meridional  plane. 

Fig.  3  displays  the  energy  spectrum  of  the  electron  differential  flux  at  L=1.2  and 
compares  it  with  the  typical  energy  spectrum  from  the  NASA  AE8MIN  empirical 
model  [Vette,  1991].  The  energy  spectrum  scales  approximately  as  E~5  with  5  ~ 
2  —  4  depending  on  the  energy  range. 

Fig.  4  shows  the  L-shcll  dependence  of  the  differential  flux  at  a  fixed  energy  of 
1.5  MeV  for  both  the  physical  model  and  the  NASA  model.  Here  the  flux  shows  a 
broad  peak  at  approximately  L=4  in  agreement  with  the  NASA  AE8MIN  empirical 
model. 

Fig.  5  gives  the  dependence  of  the  differential  flux  on  equatorial  pitch  angle  at 
L=4  and  at  a  fixed  energy  of  1.5  MeV  for  the  model  and  NASA  empirical  model. 
The  loss  cone  is  evident  at  small  pitch  angles  and  is  consistent  with  the  NASA 


empirical  model. 
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4.  Summary 

We  have  developed  a  first-principles  time-dependent  model  for  the  electron  and 
ion  radiation  belts.  The  model  consists  of  collisional  energy  loss,  collisional  pitch  an¬ 
gle  scattering,  wave-particle  induced  pitch  angle  scattering,  energy  diffusion,  cross 
energy-pitch  angle  diffusiion,  radial  L  convection,  azimuthal  (local  time)  convec¬ 
tion,  and  radial  diffusive  effects.  The  model  equations  are  in  the  form  of  a  modified 
Boltzmann  Fokker-Planck  model  and  treat  the  range  of  L-shell,  pitch  angle,  and 
energy  of  1  <  L  <  7,  3°  <  a  <  90°,  and  0.1  MeV  <  E  <  5MeV. 

The  model  has  been  applied  to  parallel  propagating  R-mode  whistler  turbulence. 
Quantities  derived  from  these  initial  studies,  i.e.,  omnidirectional  electron  flux, 
energy  spectra,  L-shell  electron  flux  dependence,  pitch  angle  diffusion  as  a  function 
of  pitch  angle  compare  favorably  with  those  derived  from  the  NASA  empirical  model 
AE8MIN. 

In  the  future  we  will  consider  other  wave-induced  pitch  angle  and  energy  diffu¬ 
sion  sources,  e.g.,  Alfven  waves  and  electromagnetic  ion  cyclotron  waves. 
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